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Abstract 

In this paper we analyze a simple mathematical model which describes 
the interaction between the proteins p53 and Mdm2. For the statio- 
nary state we discuss the local stability and the existence of the Hopf 
bifurcation. Choosing the delay as a bifurcation parameter we study 
the direction and stability of the bifurcating periodic solutions. Some 
numerical examples and the conclusions are finally made. 
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1. Introduction 

In every normal cell, there is a protective mechanism against tumoral 
degeneration. This mechanism is based on the P53 network. P53, also known 
as "the guardian of the genome", is a gene that codes a protein in order to 
regulate the cell cycle. The name is due to its molecular mass: it is a 53 kilo 
Dalton fraction of cell proteins. Mdm2 gene plays a very important role in 
P53 network. It regulates the levels of intracellular P53 protein concentration 
through a feedback loop. Under normal conditions the P53 levels are kept 
very low. When there is DNA damage the levels of P53 protein rise and if 
there is a prolonged elevation the cell shifts to apoptosis, and if there is only 
a short elevation the cell cycle is arrested and the repair process is begun. 
The first pathway protects the cell from tumoral transformation when there 
is a massive DNA damage that cannot be repaired, and the second pathway 
protects a number of important cells (neurons, myocardic cells) from death 
after DNA damage. In these cells first pathway, apoptosis, is not an option 
because they do not divide in adult life and their importance is obvious. Due 
to its major implication in cancer prevention and due to the actions described 
above, P53 has been intensively studied in the last two decades. 

During the years, several models which describe the interaction between 
P53 and Mdm2 have been studied. We mention some of them in references [2] , 
[3], [4], [6], [7], [8], [9], [10]. This paper gives a mathematical approach to the 
model described in [10]. The authors of paper [10] make a molecular energy 
calculation based on the classical force fields, and they also use chemical 
reactions constants from literature. Their results obtained by simulations 
in accordance with experimental behavior of the P53-Mdm2 complex, but 
they lack the mathematical part which we develop in this paper. We analyze 
the Hopf bifurcation with time delay as a bifurcation parameter using the 
methods from [1], [5], [7]. 

The paper is organized as follows. We present the mathematical model in 
section 2 and the existence of the stationary state is study. In section 3, we 
discuss the local stability for the stationary state of system (2) and we inves- 
tigate the existence of the Hopf bifurcation for system (2) using time delay 
as the bifurcation parameter. In section 4, the direction of Hopf bifurcation 
is analyzed according to the normal form theory and the center manifold 
theorem introduced by Hassard [5]. Numerical simulations for confirming 
the theoretical results are illustrated in section 5. Finally, some conclusions 
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are made. 



2. The mathematical model and the stationary state 



The state variables are: yi(t),y2(t) the total number of p53 molecules and 
the total number of Mdm2 proteins. 

The interaction function between P53 and MDM2 is / : IR+ — > 1R given 
by [10]: 



f{VuV2) = \{y\+V2 + k - sj{yi + y 2 + k) 2 -A yi y 2 ). (1) 



The parameters of the model are: s the production rate of p53, a the 
degradation rate of p53 (through ubiquitin pathway), and also the rate at 
which Mdm2 re-enters the loop, b the spontaneous decay rate of p53, d the 
decay rate of the protein rate Mdm2, k\ the dissociation constant of the 
complex P53-Mdm2, c the constant of proportionality of the production rate 
of Mdm2 gene with the probability that the complex P53-Mdm2 is build. 
These parameters are positive numbers. 

The mathematical model is described by the following differential system 
with time delay [10]: 

Vl (t) =s- af( yi (t),y 2 {t)) - b yi {t), ^ 
Vi{t) =cg(y 1 (t-T),y 2 (t-T)) -dy 2 (t), 

where / is given by (1) and g : IR+ — > 1R, is 

9{yi,y2) = -r— 77 v ( 3 ) 

h + yi- j{yi,y 2 ) 

For the study of the model (1) we consider the following initial values: 

yi(O) = ip 1 (8),y 2 (O) = <P2(0),6e[-T,O], 
with (pi,ip 2 : [— r, 0] — > 1R+ are differentiable functions. 

In the second equation of (2) there is delay, because the trancription and 
translation of Mdm2 last for some time after that p53 was bound to the gene. 

The stationary state (t/io, y 2 o) G IR + is given by the solution of the system 
of equations: 
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s-af(y 1 ,y 2 )-by 1 = 0, ,^ 
cg(yi,y 2 ) -dy 2 = 0. 

From (1), (3) and (4) we deduce that the stationary state can be found 
through the intersection of the curves: 

V2 = /i(yi), yi = Kyi), (5) 

where /i, f 2 : IR+ — > IR are given by: 

c(y 1 (a + b) - s) 



/i(yi) 
.M2/2) 



d(&ia + (a + b)y\ — s) ' 

(s — byi){ka + (a + &)?/i — s) 

a((o + - s) 



(6) 



Proposition 1. TTie function fi, f 2 from (6) have the following proper- 
ties: 

f s 

(i) fi is strictly increasing, f 2 is strictly decreasing on 



a + b'bj' 

(ii) There is a unique value y 10 G ( -, - ) so that /i(yio) = f 2(1/10), 

\a + b bj 

where y w is the solution of the equation (p(x) = 0, from ^ — — -, where 

f(x) = ac((a+b)x — s) 2 — d(k 1 a+ (a+b)x — s) (ka+ (a+b)x — s)(s — bx) . (7) 
Proof, (i) From (6) we have: 

(a + b)acki 



d{k ia + {a + b) yi -s) 2 ' (R) 
,,, N b bk k(a + b)(s-b yi ) W 



a (a 



+ b)y 1 - s ((a + b)yi - s) 2 ' 



For yi G ( — — -, - J, the relations (8) lead to f[(y±) > 0, f 2 (yi) < 0, so that 
\a + b bj 

fi is strictly increasing and f 2 is strictly deceasing. 

(ii) By (i) there is y w G so th at /i(t/i ) = h(yva). From (6) 

and (7) it results that: 

<p(yi) = fi2(yi)h(yi), 
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where 

h(yi) = ad((a + b)y 1 - s)((a + b)yi + ha - s) 

fn(yi) = h{yi) - hiyi). 

On ( tit)) the functions f\ 2 and h are strictly increasing and conse- 

\a + b bj 

s a 3 dkkis 

quently ip is strictly increasing. Because tpl -J = < and 

a + o a + o 

fij) = To > we can conclude that equation (p(x) = has a unique 
o b 2 

solution |/io on 



s s 



a + b 1 b 

3. The analysis of the stationary state and the existence of the 
Hopf bifurcation. 



We consider the following translation: 

y 1 = x 1 + t/io, 2/2 = x 2 + J/20 

and system (2) can be expressed as: 

Xi(t) = s - af(xi(t) + 2/io, x 2 {t) + 2/20) - b(yi(t) + y 10 ), 
x 2 {t) = cg(xi(t - t) + 2/10, x 2 (t - t) + y 20 ) - d(x 2 (t) + y 2 o)- 



(9) 



System (9) has a unique stationary state (0,0). To investigate the local 
stability of the equilibrium state we linearize system (9). Let Ui(t) and u 2 (t) 
be the linearized system variables. Then (9) is rewritten as: 



where 



U(t) =AU(t)+BU(t-r), 

/0 0\ 

c7io c7oi 



-(b + ap w ) 




B= 



(10) 



(11) 



V 
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with U(t) = (ui(t), u 2 {t)) T and p 10 , p i, 710, 7oi are the values of the first 
order derivatives for the functions: 

f( x , V) = \{x + y + k- ^{x + y + k) 2 - Axy) 
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and 

«(*.»)- x ~ f(x ' y) 



fa + x- f(x,y) 

evaluated at (2/10,2/20)- 

The characteristic equation corresponding to system (10) is A(A, r) = 
det(XI -A- e- Xr B) = which leads to: 



X 2 +p 1 X + p -(q 1 X + q )e~ XT = 0, (12) 



where 



Pi = b + d + apio, p = rf(6 + apio), ft = cyoi, 
?o = c7oi(& + apio) - ac/9 i7io- 

If there is no delay, the characteristic equation (12) becomes: 

A(A,0) = A 2 + (p 1 -g 1 )A + p -?o. (13) 
Then, the stationary state (0, 0) is locally asymptotically stable if 

Pi-qi> 0, p -qo> 0. (14) 

When r > 0, the stationary state is asymptotically stable if and only if 
all roots of equation (12) have a negative real part. We are determining the 
interval [0, r ) so that the stationary state remain asymptotically stable. 

In what follows we study the existence of the Hopf bifurcation for equation 
(10) choosing r as the bifurcation parameter. We are looking for the values 
r of t so that the stationary state (0, 0) changes from local asymptotic 
stability to instability or vice versa. We need the pure imaginary solutions of 
equation (12). Let A = ±io>o be these solutions and without loss of generality 
we assume us > 0. Replacing A = iu> and r = r in (12) we obtain: 

q cosuj T + qi uj Q sinujQTQ = p - uJq 
qosinuj T - qxUJoCOSUoTo = -u>oPi, 

which implies that 
t = — (2k + l)7i+arcsin -, = + arcsin- 



UJ 



G°o-^o) 2 +^6Ti \j(Po-ul) 2 +ulp{j 



(15) 
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where loq is a solution of the equation: 

+ (~Pi - 2po + ql)oJ 2 +p 2 - ql = 0. 



Now we have to calculate Re ( ^ ) evaluated at A = iu and r = r . We 

Ur J 



have: 



where 



and 



— U=^o,r=r = M + ?A/ 

+ 



N -qj^l + ^(gpgi - gigj +7p(2pog| - p\q\ - ql)) 

11 + 11 

^ii-piQo-^poqoqi + pfgoqi - PoPiqf+M-piqo-qiPo+ 2 Poqo)) 



+ 

with 



11 + 11 

uJoj-poPigj +plqoqi-r plql) 
H + % 

(17) 

h = -qiuj^+p^o-qxpo + r (-giPiWo -qo^i+QoPo), 
/ 2 = 2a; ? +7"o(-?i^o+Po?i^o+Pi?o^o)- 

We conclude with: 

Theorem 1. If there is no delay, under condition (14) system (10) has 
an asymtotically stable stationary state. If r > andp\q\ — q\p\ — 2p <?o > 

then there is r = r given by (15) so that Re ( — - J > and therefor 



a Hopf bifurcation occurs at (2/10,2/20)- 



/ A=io;o,t=to 



4. Direction and stability of the Hopf bifurcation 

In this section we describe the direction, stability and the period of the 
bifurcating periodic solutions of system (2). The method we use is based 
on the normal form theory and the center manifold theorem introduced by 
Hassard [5]. Taking into account the previous section, if r = tq then all roots 
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of equation (11) other than ±iuo have negative real parts, and any root of 
equation (11) of the form A(r) = a(r) ±iu>(r) satisfies a(ro) = 0, uj(tq) = Wq 

and — ; 7^ 0. For notational convenience let r = r + p, p G 1R. Then 

dr 

p — is the Hopf bifurcation value for equations (2). 

The Taylor expansion at (?/io, I/20) of the right members from (9) until the 
third order leads to: 

x(t) = Ax(t) + Bx(t - t) + F(x(t),x(t - r)), (18) 

where x(t) = (xi(t),x 2 (t)) T , A, B are given by (11) and 

F(x(t),x(t - t)) = (F 1 (x(t)), F\x{t - r))) T , (19) 

where 

F 1 (x 1 (t),x 2 (t)) = ~[p2oxj(t) + 2p 11 x 1 (t)x 2 (t) + p 2X 2 2 (t)}- 

- I [p 30 xl (t) + 3p 21 xl (t)x 2 (t) + 3pi 2 xi (t)a| (i) + p 03 x^ (*)] 

F 2 (a:i(t-r),x 2 (t-r)) = ^ 20 xl(t - t) + 2 lllXl (t - r)x 2 (t - t)+ 

+ 702X2 (*- r)]+ 

+ ^ [730^1 (< + 3721X1 (t - r)x 2 (t - r) + 
o 

+ 3712X1 (t - r)xl (t - r) + 703^1 (t - r)\ 

and pjj, 7jj, i,j = 0,1,2,3 are the values of the second and third order 
derivatives for the functions f(x,y) and g(x,y). 

Define the space of continuous real- valued functions as C = C([— r , 0], 1R 4 ). 

In r = t + p, p G H, we regard p as the bifurcation parameter. For 
$ G C we define a linear operator: 

LQu)$ = A$(0) + B$(-r) 

where A and B are given by (11) and a nonlinear operator -F(/i, $) = F($(0), 
$(-t)), where F($(0), $(-r)) is given by (19). According to the Riesz 
representation theorem, there is a matrix whose components are bounded 
variation functions, r](6, p) with 9 G [—To, 0] so that: 



L(/x)$= / d v (e,p)<f>(9), [-7-0,0]. 

-TO 
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For $ G (^([-TcO^IR 4 ) we define: 



dO ' 



9e[-T ,0) 



J dri(t,ri<!>(t), 6 = 0, 

-TO 

0, 9e[-T ,o) 
F(/i,$), 6> = 0. 
We can rewrite (18) in the following vector form: 

tit = A(n)ut + R(fj)ut (20) 

where u = (ui, u 2 ) T , u t = u(t + 6) for 6 G [—To, 0]. 

For * G ^([0, r ],IR* 4 ), we define the adjunct operator A* of .4. by: 



A*V(s) = < 



s G (0, r ] 



/ dr) T (t,0)i>(-t), s = 0. 



V -r 

We define the following bilinear form: 

f-0 



< m >= ^ T (o)$(o) - / / ^ T (e - e)dri(o)<KS)dZ, 

J -TO J£=0 

where r)(6) = r}(8,0). 

We assume that ±icu are eigenvalues of .4(0). Thus, they are also eigen- 
values of A*. We can easily obtain: 

Q>(9)=ve Xie , 0e [-7b,0] (21) 

where v — (v±, v 2 ) T , 

vi = -ap 0l ,v 2 = Ai + b + apio, 
is the eigenvector of .4(0) corresponding to Ai = iujq and 

m(s) = we X2S , sG[0,oo) 
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where w = (wi,W2), 

w 2 cyi e~ XlT 1 

Wl = 7 , x ~ ,W 2 = ~, 

b + Ai + ap w rj 

„ = ^f"" + v, - ™*+<™'» { - Ta X ie -^ - 1 + )] 
o + A 2 + apio Af 

is the eigenvector of .4.(0) corresponding to A 2 = — icJo- 
We can verify that: < #(s),$(s) >= 1, < *(s),l>(s) >=< ^(s), $(s) >= 

0, < tf(s),$(s) >= 1. 

Using the approach of Hassard [5], we next compute the coordinates to 

describe the center manifold Q at fi = 0. Let u t — u t (t + 9), 9 e [—To, 0) be 

the solution of equation (20) when fi — and 

=< *, ^ >, «;(*, 0) = u t {9) - 2Re{z(t)<S>(6)}. 
On the center manifold Q , we have: 

w(t,9) = w(z(t),z(t),9) 

where 

^2 ^2 ^3 

iu(z, 5, 0) = w 2 q{9)— + iu n (0)zz + wq 2 {9)— + w 30 (9)— + ... 

in which z and z are local coordinates for the center manifold Qq in the 
direction of \& and \l/ and w 02 (9) = w 20 (9). Note that w and M t are real. 
For solution u t £ f2 of equation (20), as long as /i — 0, we have: 

i(t) = Aiz(t)+0(z,;z) (22) 

where 

z) = *(0)F{w(z(t), z(t), 0) + 2ite(z(f)$(0))) = 

= ^20^ + 9liz(t)z(t) + 302^" + g21 W 2 + • • • 

where 

#20 = F 2 \w x + F 2 2 w 2 , c/ii = F^wi + F^w 2 , g 02 = F^w-y + F$ 2 w 2 , (23) 
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with 



F 2 o = -a(p 2 ovf + 2p 11 v 1 v 2 + p 2vj, 
Fi = c^vfe- 2 ^ + 7o2V 2 2 e- 2Air + 2 lllVl v 2 e- 2X ^], 
F ii = ~a[P2oViVi + pii(^iv 2 + viv 2 ) + P02V2V2], 
F11 = c[72o^i% + 711(^1^2 + viv 2 ) + 702U2V2], 



pi _ Z?l Z?2 _ p2 
^02 ~~ ^20' ^02 ~~ ^20' 



and 
where 



021 = ^21^1 + F 2 2 1^2 (24) 



21 



-a[p 20 (2v 1 «; 1 1 1 (0) + ui«4(0)) + 2 Pll (*Wi(0)+ 

+ + ^ l(0)) + ^»( 2 ^i(°) + ^<(0))- 

P30^i + 2p2i^i^2^i + 2p 12 v 1 v 2 v 2 + Po3viv 2 + p 2 \v\v 2 + p\ 2 v x v 2 2 \ 

F 21 = c[^ 20 (2v 1 w 1 u (-r)e- Air + v 1 «; 2 1 (-r)e A ^) + 2 lll (v l w\ 1 (-T)e-^ T + 

+ VlW ^ )eXl \ V ^^ 

+v 2 w 2 (-r)e XlT ) + 7 30 ^i 2 Vie" AlT + 72i(2viVit> 2 e- Air + vjv 2 e- XlT ) + 
+~f 12 {2v 1 v 2 v 2 e~ XlT + Vivle~ XlT ) + 7o3^2e~ AlT ]. 

The vectors w 20 (9), wu(9) with 6 e [— r, 0] are given by: 



(25) 



/ /i\ #20 Ai6» #02 - A 2 9 TTi 2Ai0 

u»2o(p) = — —we 1 —ve 2 + E ± e 1 

Ai 0A1 

//n 9ll Ai0 #H - Ao6» 1 771 

= —we 1 — — ve 2 + £/ 2 
Ai Ai 

where 

E l = -(A + e~ 2X ™B - 2X 1 I)- 1 F 20 , E 2 = -{A + B)~ 1 F 11 , 

where F 20 = (F 2 \,F 2 ) T , F u = (F^F^f. 

Based on the above analysis and calculation, we can see that each in 
(23), (24) are determined by the parameters and delay of system (2). Thus, 
we can explicitly compute the following quantities: 

Ci(0) = ^(010011 " 2 l^i| 2 " 5 W 2 ) + 9 -f, 

RejCM) rp ImidjO)) + p 2 N 9n , r([) ^ (26) 
P2 = 77 , T 2 = ,p 2 = 2Re{C 1 {0)), 

IVl LOq 
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where M and N are given by (16) and (17). 
In short, this leads to the following result: 

Theorem 3. In formulas (26), /i 2 determines the direction of the Hopf 
bifurcation: if fi 2 > 0(< 0), then the Hopf bifurcation is supercritical (sub- 
critical) and the bifurcating periodic solutions exist for r > t (< t ); (3 2 
determines the stability of the bifurcating periodic solutions: the solutions 
are orbitally stable (unstable) if f3 2 < 0(> 0); and T 2 determines the pe- 
riod of the bifurcating periodic solutions: the period increases (decreases) if 
T 2 > 0(< 0). 



4. Numerical example. 

In this section we find the waveform plots through the formula: 

X(£+0) = z(£)$(0)+z(t)i.(0)Hi^ 

where z(t) is the solution of (20), <&(#) is given by (21), w 2 o(9), wn(6) , wq 2 (9) = 
w 20 (9) are given by (25) and X = (yio,y 2 o) T is the equilibrium state. 

For the numerical simulations we use Maple 9.5. and the data from [10]: 
the degradation of p53 through ubiquintin pathway a = 3 x 10 _2 sec~l, the 
spontaneous degradation of P53 is b = 10~ 4 sec _1 , the dissociation constant 
between P53 and Mdm2 gene is k\ = 28, the degradation rate of Mdm2 
protein is d = 10~ 2 sec _1 , the P53 protein production rate is S — O.Olsec -1 
and the production rate of Mdm2 is c = lsec -1 . For this date we consider 
the different values for the constant k. By changing instantaneously the 
dissociation constant k the response of the system is different. 

For k = 17.5 we obtain: the equilibrium point yw = 1.674122637, y 2 o = 
4.587857801, the coefficients which describe the limit cycle: /i 2 = 0.002340098338, 
(3 2 = -0.000001452701020, T 2 = 0.00007793850775 and u = 0.04577286901, 
r = 60.52296388. Then the Hopf bifurcation is supercritical, the solutions 
are orbitally stable and the period of the solution is increasing. The wave 
plots are displayed in figl and fig2: 
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Fig.l. (t,yi(t)) 



Fig.2. (t,y 2 (t)) 
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For k = 120 we obtain: the equilibrium point yio = 3.853801769, y 2 o = 
11.20502079, the coefficients which describe the limit cycle: /i 2 = -8. 219839362- 
10- 7 , / 9 2 = 2.0434666644-10- 10 , T 2 = 2.243699328- 10" 7 and oo = 0.02923152517, 
ro = 92.68683554. Then the Hopf bifurcation is subcritical, the solutions are 
orbitally unstable and the period of the solution is increasing. The wave 
plots are presented in fig3 and fig4: 



Fig.3. (t,yi(t)) 



Fig.4. (t,y 2 (t)) 



l\ h (\ 
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For k = 1750 we obtain: the equilibrium point y±o = 14.88816840, 
1/20 = 34.27918463, the coefficients which describe the limit cycle: \i 2 = 
1.104273140 • 10 -11 , 02 = -6.4139104 • 10" 16 , T 2 = 1.924908056 • 10" 11 and 
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lu = 0.01423761906, r = 174.4149631. Then the Hopf bifurcation is supercri- 
tical, the solutions are orbitally stable and the period of the solution is in- 
creasing. The wave plots are given in fig5 and fig6: 




5. Conclusions. 

For the present model, we obtain an oscillatory behavior similar with the 
findings in [10], according with the qualitative study. 

We have proved that a limit cycle exists and it is characterized by the 
coefficients from (26). For different values for the equilibrium constant k, 
in Section 4 we obtain stable or unstable periodic solutions with increasing 
periods, via a Hopf bifurcation. 
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